Reply to Comment on: "Are stress-free membranes really 'tensionless' ? 



Friederike Schmid 

Institute of Physics, JGU Mainz, D-55099 Mainz, Germany 



PACS numbers: 87.16.dj, 68.03.Kn, 68.35.Md 



(N ■ 

o : 



o 



-a 
c 

o 
o 



> 
o 

O 

(N 



X 



Fournier and Barbctta state the central message in 
their paper as follows (see conclusion of Ref. Q]): " We 
showed that r differs from the tensionlike coefficient r 
of the fluctuation spectrum and we unveiled the correct 
way to derive r from the free energy." Here r is the 
frame tension, called a ( in my work, and r is the fluctu- 
ation tension, called a nuo in my work. In Ref. [H, I argue 
that the arguments of Fournier and Barbetta - as well 
as those of other previous authors who came to the same 
conclusion 0, 0] - are inherently inconsistent, due to the 
fact that they are based on a theory which has been lin- 
earized with respect to a small parameter (A — A p )/A p , 
yet predict an effect which is nonlinear in this parame- 
ter. (Here A is the membrane area and A p the projected 
membrane area.) Fournier now claims that the argu- 
ments in Ref. [l| is really based on an expansion in k B T, 
which would be the usual expansion in a diagrammatic 
field-theoretic treatment of the problem. 

Let us therefore consider the expansion in k B T. The 
starting point is the Helfrich Hamiltonian Whsis** hi 
Monge gauge, i.e., a gauge where the coordinates (x,y) 
are defined by projection onto the plane of the frame and 
the membrane is parametrized by the normal distance 
h(x, y) onto this plane. In the field-theoretic treatment, 
% H oifrich i s expanded about a planar reference state, yield- 
ing i 
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Here o~ is the bare tension, which coincides with the 
internal tension a int in a system with fixed number of 
lipids, and k is the bending rigidity. The Hamiltonian 
Ho is Gaussian and %' subsumes the nonlinear terms. 
The k B T expansion is based on the full nonlinear Hamil- 
tonian %, where H' is treated as a perturbation about %q 
in a diagrammatic scheme. Further corrections come in 
through the nonlinearity of the measure T>[h] Q. Within 
this approach, some quantities can be evaluated up to 



first order k B T without having to actually consider non- 
linear corrections to "Ho- The frame tension a t is probably 
such a quantity, and therefore, Eq. (4) in Ref. |2j indeed 
gives the correct leading correction for a { /a int in a k B T 
expansion. The fluctuation tension <7 fluc , however, is not 
such a quantity and the calculation of the k B T or- 
der involves the calculation of first-order loop diagrams. 
Such calculations can be very tricky [(| and have not 
been attempted in Rcfs. [l|, 0, If Fournier and Bar- 
betta meant to show o> ^ er fluc by an expansion in powers 
of k B T, then their calculation not just inconsistent, it is 
incomplete. In that case, they should have finished the 
one- loop calculation for er fluc before makingany claims. 

We conclude that the arguments of Refs. [lj, |j clearly 
cannot be justified by an expansion in k B T or the cor- 
responding dimensionless quantity e = k B T/n. In- 
stead, the authors Refs. El, 0, have simply replaced 
% = %q, which implies (among other) omitting higher 
order terms in (V/i) 2 -C 1 and settina <7 nilc = cr int . This is 
the approximation examined in Ref. y. From the relation 
— (V/i) 2 dA p , one gets locally 
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with rj(x,y) = d(A — A p )/dA p , hence (Vft.) 2 <C 1 implies 
i] <C 1. The global average of r] is fj = (A — A p )/A p , 
which is thus a small parameter in this approximation: It 
neglects terms that are not linear in fj. The "expansion" 
is not systematic, because other terms (higher orders of 
higher derivatives of h) are neglected as well, but this is 
not important for our argument. 

It is important ot note that the parameters fj = 
{A— A p ) I A p and e = k B T / k can be varied independently. 
This is physically feasible, since A p can be controlled ei- 
ther directly or by tuning the frame tension, independent 
of the temperature k B T. The approximations fj <C 1 and 
e <C 1 arc thus not equivalent. On the one hand, capil- 
lary wave Hamiltonians || that ignore bending terms - 
which corresponds to setting k— >0ore— >oo - have been 
extremely successful in describing the properties of liq- 
uid/liquid interfaces at large wavelengths. On the other 
hand, membranes with fixed number of lipids and ap- 
proximately fixed area per lipid can be studied at fixed 
projected area. This is actually a common setting in sim- 
ulations. 
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By appropriate Legendre transforms, one can calcu- 
late the free energy of the Gaussian model Hq in such a 
(JV, A, A p ) ensemble: 



F(N,A,A P ) 
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(4) 



The frame tension and the fluctuation tension tension 
can be calculated via o~ { = dF/dA p and a nnc = cr int = 
—dF/dA, and the results are of course the same as those 
presented in Refs.[l|-[i]for other ensembles. Nevertheless, 
the results from the Gaussian model clearly cannot be 
trusted at order ((A — Ap)/A p ) 2 . Nonlinear effects will 
become important even at small temperatures if (A — 
Ap)/A p is large. 

I wish to stress once more that this whole controversy is 
not about the relation between the frame tension and the 
internal tension, but about the fluctuation tension. One 
should give Farago and Pincus credit for having been the 
first to derive the relation (4) in Ref. [H between cr f and <7 int 
for compressible membranes with fixed number of lipids 
Q. As Fournier correctly points out, this result gives 
most likely the correct leading order in a diagrammatic 
expansion in powers of k B T/n. However, <x int is a rather 
uninteresting quantity, since it can neither be controlled 
nor measured. Farago and Pincus recognized that their 
result does not carry over to the fluctuation tension, and 
gave a very general argument why the fluctuation ten- 
sion should equal the frame tension 1(| U | , which solely 
relies on the requirement of " rotational invariance", i.e., 
gauge invariance. Their reasoning is similar to a classic 



argument by Cai et aZ.[6j], who showed o~ aac = a t for in- 
compressible membranes with variable number of lipids. 
The conclusion that gauge invariance leads to cr tluc = <j { 
has ver y re cently been corroborated by numerical simu- 



lations 11 



Notwithstanding, the highly accurate simulations of 
Ref.0 suggest that er fluc should be slightly renormalized, 
Cauc = o~ { {Ap/A). This result is in line with model-free 
thermodynamic considerations on the relation between 
different tension parameters in vesicles [12]. Whether 
and how it can be reconciled with the general arguments 
for <7 fluc = er t quoted above still remains to be elucidated. 
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